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Abstract The experimental muon source on China Spallation Neutron Source (CSNS) is expected to be a high 


intensity (10° y'/s) surface muon source with a small beam spot of 4-cm diameter. For a practical application of this 


muon source, we are devoting to develop the first pulsed uSR spectrometer in China. In this paper, the performance of 


plastic scintillators in the uSR spectrometer is studied by Monte Carlo simulation. The processes such as positron 


energy deposition, scintillation photons production, light propagation and photon-electron conversion are carefully 


considered. According to the results, an optimal dimension of the plastic scintillator is proposed using for our future 


spectrometer, which has a long-strip shape with the dimension variation range of 50-60 mm length, 5-8 mm height, 


and 10-12 mm width. Finally, we can build a spectrometer with a count rate up to 10fe*/s by 100-120 forward and 


backward segmental detectors in total. The simulation could serve as an important guide for spectrometer construction. 
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1 Introduction 


China Spallation Neutron Source (CSNS) is an 


accelerator-based facility, including a proton 
accelerator, a neutron target station and multiplex 
neutron scatter spectrometers. An intensity neutron 
source will be produced by a 1.6-GeV proton beam 
bombarding a high Z metal target!'7! In the phase I, 
120-kW proton beam, neutron beam intensity of 
2.0x10'° n/em’-s, and 25-Hz synchronic frequency 
will be achieved. Besides serving as a neutron source, 
CSNS can also be used as a potential muon sourceľ!, 
The Experimental Muon Source (EMus)**), which is 
proposed for Phase I, is expected to perform a surface 
muon source with beam intensity of ~10°u"/s, beam 
spot of 4-cm diameter, and 1-Hz pulse frequency. 
Using this new muon source, a spectrometer for uSR 
spectroscopy (Muon Spin Rotation, Relaxation, 
Resonance, etc)!*”! is essential to be developed now. 


Fig.l schematically shows a pulsed uSR spectrometer. 
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Nearly 100% polarized surface muons are injected into 
specimen and decay into positrons asymmetrically 
distributed in the full 42 solid angle. A muon counter 
before the uSR spectrometer is used to count incoming 
muons as a Start-signal, and two sets of segmental 
detectors located forward and backward to the sample 
are used to detect the positrons as stop-signals. A 
lifetime of muons in materials could be measured by 
registering the start- and stop-signals. Finally, we 
could obtain the positron spatial distribution as time 
goes to explore the internal structure and dynamics in 


condensed matters under various extreme conditions. 


Forward and 


Fig.1 Schematic of the uSR detection system. Positrons from 
muon decay are asymmetrical distributed into space with an 
angular distribution: P(@)=1+Acos(a). 
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Type of uSR spectrometer is dependent on the 
time structure of muon source: continuous or pulsed"), 
Benefiting from the development of detection and 
superconductor techniques, people are committing 
themselves to constructing 


some novel wuSR 


spectrometers with improved performance. The 
photodiode detectors with small dimensions, high 
detection efficiency, and less magnetic sensitiveness 
are used to replace the traditional PMT in order to get 
higher resolution of 100 ps7). Superconductors are 
used to supply a high external magnetic field up to 5 T 
for the special uSR study!” Silicon microstrips are 
served as position-sensitive detectors for the study of 
spatial-temporal uSR spectroscopy"), New data 
acquisition systems are developed for vast scintillators 
detector system!*), The realms of uSR spectroscopy 
will be dramatically broadened by these advanced 
pulsed wuSR 
construction, we have to utilize many small detectors 


techniques. For a spectrometer 
in order to detect amounts of positrons in a short time 
interval of around 2.2 us. The sizes of scintillator 
detector play a significant role in structure designing, 
optical coupling, light transmission and signal 
processing etc. Unfortunately, the performance of the 
positron detectors are impossibly conducted using a 
variable energy positron beam of 0-53 MeV, one only 
decide the scintillator parameters by experience. So, 
we performed some Monte Carlo simulation to study 
the plastic scintillator in the uSR detector system, and 
other simulations for uSR spectrometer construction 
could be found in Refs. [13,14]. 

In Section 2, the simulation configuration, 
parameters and approaches are briefly introduced. 
Many physical quantities are studied, such as the 
energy deposition and the muons energy loss rate (K) 
in a scintillator, the transition probability (P) of 
positron energy to scintillation photons, light yield (R) 
and the light transmission efficiency (7) in some 
mediums. The simulation results were analyzed and 
discussed in Section 4. Finally, a preliminary proposal 
for our future spectrometer construction is given based 


on our simulation. 


2 Simulation, configuration and methods 


The physical process is to be simulated as follows: the 
high intensity surface muons up to 10°u'/s are injected 


into specimen, stopped and decayed into positrons, the 
positrons are detected by transporting and collecting 
the scintillation photons produced by the ionization 
loss of energy in scintillators. In this process, the 
electromagnetic interactions among particles and 
materials and optical processes can be simulated by a 
Monte Carlo code — Geant4!'*'®!, This tool allows us 
to set many flexible and modular simulation 
conditions, such as the detector dimensions, the 
reflector configuration used for scintillator wrapping, 
the properties of optical mediums with an optical 
coupling and the realistic scenario of PMT response. 
Fig.2(a) shows a single channel for positron detection 
by a plastic scintillator in the uSR spectrometer system. 
Simulation was carried out using the long-strip EJ-200 
plastic scintillators with dimensions of 30-100 mm 
length, 5-15 mm height and 5-15 mm width. The 
main characteristics of the scintillator are summarized 
in Table 1(a). Length (L) and width (W) will determine 
the count rate of the single-channel detector or the 
integral spectrometer. Meanwhile, light transmission 
efficiency and positron energy deposition are evidently 
dependent on length (L) and height (H). In our 
simulation, the photon exit surface of the scintillator 
HxW is connected with a light guide end of 5 mmx5 
mm using the silicone optical grease of 1.5-refractivity, 
while the other surfaces are covered with a Teflon 
reflective coating. A 1-meter-long clear fiber BCF-98 
served as light guide is used to connect the scintillator 
to transport the scintillation light to a PMT (ET 9813B 
(Table 1b). Detection and measurement can be free 
from the high magnetic filed of several thousands 
Gauss in the central area of spectrometer. This clear 
light guide, which has double claddings of 0.5 mm 
thick Acrylic inner and 0.1 mm thick Fluor-Acrylic 
outer, is reported in Table l(c). The maximum 
incidence angle of the total- reflection light in the fiber, 
which can be defined as @max with respect to the long 
axis of the fiber, is given by singmaxV[(n1-—n3 Vno], 
where, nı and n; are the refractive of the core material 
and the outer cladding of the fiber. nọ is the refractive 
of the silicon grease coupled with the fiber ends. The 
max around 30 degree calculated by the parameter 
values can be obtained. Therefore, a fraction of 
photons exiting out the scintillator will be inevitably 
lost at the front end of light guide. The last part of our 
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simulation in Fig.2(a) is the 


conversion by modeling the PMT cathode realistically, 


photon-electron window connected with the exit end of light guide and 
blue-green sensitive bialkai photocathode. 


including a 5 mm thick, 51 mm diameter quartz 


(a) (b) 


PMT Window 
and Cathode 


/ Light Guide 


Plastic Scintillator 


Surface muon 


Backward Segmental Detectors 


NS Forward Segmental 
Detectors 


Fig.2 (a) Sketch to illustrate the principle used for positron detection by plastic scintillator. Positrons from muon decay inject into 
the scintillator strip at an incident angle 0 with respect to the normal of the incident surface. (b) 3D view of segmental detectors in 


uSR spectrometer. 


Table 1 Tables with detector properties used in simulation 


(a) EJ-200 scintillator 


Base 

Atomic ratio (H:C) 

Density (g cm *) 

Refractive index 

Light output (% anthracene) 
Scintillation yield (photons/MeV) 
Emission peak (nm) 

Decay time (ns) 

(b) ET 9813B 

Refractive index 

Effective photocathode diameter(mm) 
QE peak (%) 

Gain 

Rise time (ns) 

Transit time(ns) 


(c) BCF-98 

Material 
Core Polystyrene n\=1.60 
Cladding 1 Acrylic no=1.49 
Cladding 2 Fluor-acrylic n3=1.42 


Refractive index 


Polyvinyltoluene 
11:10 

1.023 

1.58 

64 

10,000 

425 

2.1 


Numerical aperture 
/ 

0.58 

0.74 


back painted (Sigma Alpha=0.2)!"*! and polished for 


For realistically treating the interactions of the 
scintillation photons with dielectric surfaces, two 
optical reflection models—the glisur model and the 
unified model!'” were carefully set in our simulation. 
The surface type was dielectric_dielectric photodiodes 
for the Teflon reflector, where the unified model was 
used with the reflection type set to 100% specular lobe 
distribution. For all grease-coupled surfaces, 
dielectric_dielectric was set as surface type with the 
glisur model so as to use the Fresnel formulas for 


optics calculation. The surface finish was set to ground 


the unified and the glisur model. To keep the 
scintillator and the light guide away from the outside 
light exposure, the dielectric_metal with the glisur 
model (reflectivity: 0, quantum efficiency: 100%) is 
defined as a photon absorber attached to the outermost 
layer. In order to directly count the photoelectrons 
converted by the coming scintillation photons, 
dielectric_metal with the glisur model (reflectivity: 0.2, 
quantum efficiency: 30%) is defined as a PMT cathode 
attached behind the quartz window. 
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To estimate the positron count rate of a single 
detector or the spectrometer, we studied a traditional 
spectrometer including two symmetrical detector 
segmentations forward and backward (Fig.2b). The 
inner end surface (Hx W) of either segmental detectors 
ring are 7.5 cm away from the sample center, and the 
detector unit is 10 cm far from the ring center so as to 
save moderate space to match a compromise for 
enough sample chamber space and sufficient detection 
solid angle. 


3 Results 


Firstly, the 
scintillators with variable dimensions was studied by 


positron deposited energy in the 
perpendicularly injecting positron into the down 
surface and counting the exit positrons from the upper 
surface. The deposited energy can be calculated by 


Energy Deposit / 
N WW 
uw oOo v 


H N 
“n o 
+ + 


0 5 10 15 20 25 30 35 40 45 50 55 
Positron Incident Energy/ MeV 


subtracting the exit energy from the coming energy. 
The scintillators in this simulation have a fixed length 
of 40 mm and width of 10 mm, but variable height of 
5-16 mm. Fig.3(a) shows the deposited energy as a 
function of positron incident energy of 0-53 MeV for 
different scintillator heights. Fig.3(b) shows the energy 
loss ratio as a function of positron incident energy for 
variable scintillator heights. An empirical equation for 
energy loss ratio (K) depended on incident energy (£) 
and penetration height (H) can be obtained. 


K(E,H)=(0.12e"'°?+0.017)+(0.084e “"+0.0072)H (1) 


Where £ is the range of 5 to 53 MeV, H is penetration 
distance (mm). The constant coefficients for Eq.(1) are 
calculated by linear fitting (R°>0.98) and two 
single-exponential curves fitting (Chi*2/DoF=0) based 
on the data (Fig.3b). 


T T T T T T T T T 


o7t $ (b) | 
osl į 
0.5 
0.4} 
0.3} 
0.21 


Energy loss ratio / K 


0.1; 


10 15 20 25 30 35 40 45 50 55 
Positron Incident Energy/ MeV 


0 5 


Fig.3 (a) Deposited energy as a function of positron incident energy for nine different height scintillators; (b) Energy loss ratio as a 


function of positron incident energy for variable scintillator height. 


Figure 4 shows that the light yield of positrons 
in the scintillator was calculated with the same sizes as 
the previous deposited energy simulation at variable 
height of 5-16 mm. The vertical axis means the total 
light yield of one parametered scintillator as a function 
of positron energy of 0-53 MeV. The filled dots (a, b, 
c and d) mean the maximum light yield at some certain 
energy for different scintillator heights. Comparatively, 
the dot O is marked as the light yield of 0.511 MeV 
y-rays in the scintillator. As mentioned in Section 2, 
the scintillation photons are transported far from the 
sample room via a long light guide, avoiding the high 
magnetic field influence. Consequently, a considerable 
light loss dependent on the scintillator sizes and light 
guide is expected at the beginning of light guide. 


=. e n N W 
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Photon production per positron 


10 15 20 25 30 35 40 45 50 55 
Positron Incident Energy / MeV 


=“ 0 5 


Fig.4 Scintillation photon yield for variable scintillator under 
different incident energy. The red filled dots a, b, c and d 
indicate the approximate maximum photon production for the 
variable scintillator heights. Dot “O” means the photon 
production for the 0.511-MeV gamma rays. 
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In our case, the loss is up to 40% when setting 
an exit scintillator surface of 5 mmx10 mm (#xW) 
and an entry light guide surface of 5 mmx5 mm. For 
realistically studying, the positrons are to be injected 
into scintillators with a variable incidence angle (0) 
(Fig.2a). On the other hand, the length of the 
scintillator strip determines the @ range in our 
simulation. For example, a 50-mm-long scintillator 
strip means that 6 range of 38 degree to 50 degree was 
calculated according to the spectrometer configuration. 
Fig.5 shows the photoelectron count as a function of 0 
incident angle for six different lengths. The solid lines 
are linear fitting for these data. 


Photoelectron Number 


4 Se eS eS S — es — 
36 38 40 42 44 46 48 50 52 54 56 58 
Incident Angle / Degree 


300 Jad 


Fig.5 Photoelectron yield as a function of e* incident angle 
for different scintillator lengths. The lines indicate the linear fit 
for the different lengths. The correlation coefficient R? in the 
legend represents the goodness of linear fit, which is in range 
from 0 to 1. 

The expected pulsed width of surface muon 
beam at EMuS is about 100 ns smaller than the 2.2-us 
muon lifetime. We have to detect the majority of 
positrons from muon decay in 2.2 us after the muon 
beam reaching the specimen. When defining the 
Minimum Detection Time per Event (MDTE), the 
detection time for a single event in a detector unit 
should not be less than the ratio of one 2.2-us muon 
lifetime to the total counts of positrons hitting on a 
single scintillator strip. A small MDTE is important for 
us to avoid data piling up. Fig.6 shows that the MDTE 
for the different length and width of scintillators in 
simulation had an intensive surface muon beam of 
10°u"/s. Selecting the crystals with different width (8, 
10, 12, 20 mm) represents the detector ring with 
different segmental counters to be configured. For 


example, an 8-mm-width plastic scintillator used as 
detector unit means 78x2 segmentation counters 
capably integrated in our spectrometer. The other 
width size can be detailed in the inset. Additionally, 
the scintillator length is from 30 mm to 70 mm. 


— 


e.. 


a wiatnj/mm & | J 

45 Segment detector 4 
—e-—8 & 78x2 4 

40 Š —e—10 & 63x2 4 


a -4—12 & 52x2 | ] 
35 = —*—20 & 31x2 | 1 


ns 


MinmumTetection Time per Event / 


40 45 50 55 60 65 
Length / mm 


25 30 35 70 75 


Fig.6 Minimum detection time per event for different 
dimensions of scintillators. The dashed box means that the 
parameter option in this area meets the requirements of our 
optimal design. 


4 Discussion 


For a given height scintillator, we can find that 
deposited energy increases with the incidence energy 
of positrons (Fig.3a), whereas the energy loss ratio 
K(E) indicates an exponentional reduction (Fig.3b). 
But in the case of injecting a given energy positron, 
the longer path positron experience (H) is, the more 
energy deposited in scintillator is. And the function 
K(H) approximately behaves as a linear increase. 
Additionally, it is not a linear relation between the 
light yield and the incidence energy (Fig.4). For the 
scintillator height of 5 mm to 16 mm, we define the 
positrons with energy of less than 5 MeV as Low 
Energy Positron (LEP). The LEP can be completely 
stopped in the scintillator, resulting in a light yield in 
proportion to incident energy. Further, the stopped 
positrons will be annihilated with electrons, yielding 
two 0.511 
scintillator by these gamma rays will add light yield to 
our calculation. The filled dots (a, b, c, and d in Fig.4) 
give an evidence that the 0.511 MeV gamma rays do 


MeV gamma rays. Illumining the 


produce many additional photons. However, for 
positron energy of larger than 5 MeV so-called 
Medium Energy Positron (MEP), the light yield is not 
linearly increasing with incident energy. This is mainly 


due to the limited scintillation light yield in a fixed- 
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length scintillator for the Minimum Ionizing Particles 
(MIP). Consequently, its light yield R in a given 
scintillator is only dependent on the height size (H 


pR 


where, £ is in the range of 5 to 53 MeV. 

It is concluded that the decreasing energy loss 
ratio and conversion efficiency, which are opposite to 
the increasing deposited energy, exactly ensure the 
output electronic signals with close amplitudes for 
different position incidence energy. With this result, 
we can easily distinguish the MEP and background 
particles in detection, and realize positrons detection 
with nearly full energy spectrum and indiscrimination 
for energy range of 5 to 53 MeV. Avoiding the 
undesired LEP positron detection, a thin aluminium 
layer is usually placed before the scintillator!" ”, 
According to the light yield in scintillators with height 
5 mm, the light yield count is large enough up to 10°, 
for us to distinguish the background of a 0.511-MeV 
light yield (for comparison in Fig.4 filled dot O). In 
consideration of scintillator-guide coupling and the 
cost of light guide, we selected the 5-mm height 
scintillator for our simulation and future experiments. 
The distance (D) between the incident position of 
positrons and the light-exit surface of the scintillator 
significantly determinates the light transmission 
efficiency. An optimal length of the scintillator strip 
has to be found for our spectrometer design. Generally 
speaking, there is a liner relation between photons loss 
in the strip and the D. Especially for strip of larger 
than 60-mm length, the values of R? are close to 1 
(Fig.5), meaning that the photon loss in scintillator is 
well linear with the D. The large length of the 
scintillator strip is the main reason of the photon loss. 
For 80-mm length scintillator, the efficiency of photon 
collection for the positrons with @=38° (positrons 
hitting on the right end of scintillator in Fig.2a) is only 
~60% for 0= 56°(light-exit end). In this case, the final 
electronic signal amplitudes will be so evidently 
different that it is difficult for discrimination threshold 
setting. However, for the length 30 mm and 40 mm, it 
is found that the photon collections for different 
incident angles are not marked different but fluctuant. 


(—442+1880H)hv 
EK CO + 0.017) + (0.08467? + 0.0072)H |E 


mm): R(H) = —442 + 1880H. Finally, the effective 
efficiency of the energy-light conversion can be 
obtained by Eq.(2). 


(2) 


There is no clear regularity for the photon counts at 
different incident angles. The very small R (< 0.2) 
clarifies the slight relation between length and photon 
loss. The fluctuation of photon count may be mainly 
due to the different 
fluctuation or some 


incident angle, statistical 


optical physics processes. 
Consequently, the scintillator strip with length 30—60 
mm is suitable for us to control the signal difference of 
less than 20%. In this case, we can obtain an average 
photoelectron count up to 500 electrons per event, 
resulting in a voltage output up to ~1 V after the 
photoelectrons multiplier by PMT (gain 10°, load 
resistance 50 Q). However, a determinate length 
scintillator must be associated with the detection rate 
which is strongly dependent on length and width. 
From Fig.6, the selecting a scintillator strip with larger 
width and length requires a higher detection rate of 
single detection channel or integral spectrometer. The 
detection time per event can be achieved in a low 
value about 20-30 ns for the traditional detection 
equipments, including the fall time in scintillator, 
transmission time, the rising time in PMT and the 
pulsed width necessary for electronic system in total. 
With a comprehensive consideration, the scintillator 
with 50-60 mm length and 10-12 mm width is 
optimal for the positron detection in our wuSR 
spectrometer design (Black dashed box in Fig.6 shows 
this optimal result). Finally, we can obtain a desired 
performance of our future spectrometer based on the 
design in this paper: The 100-120 detection channels 
can be sufficiently segmented on the forward and 
backward detector rings, the detection rate can be up 
to 70-90 e*/s and ~1.0x10* e*/s for a single detection 
unit and whole 


spectrometer, respectively. The 


performances meet the requirements of the uSR. 


5 Conclusion 


We have used Geant4 code to study the 


electromagnetic processes and the optical processes in 
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muon lifetime detection, and investigated the suitable 
setup and the expected performance of the uSR 
CSNS EMuS. The 
parameters as deposited energy, light yield, energy- 


spectrometer at significant 
light conversion efficiency and transport efficiency 
were calculated and studied, supporting that the 
scintillator sizes play an important role in uSR 
spectrometer construction. A plastic scintillator strip 
with 50-60 mm length, 10-12 mm width, and 5 mm 
height was testified as an optimal positron detector for 
our uSR spectrometer. The high light yield and the 
light transmission efficiency in this scintillator meet 
the requirements for signal discrimination and uSR 
a wuSR 
detection 


spectroscopy measurement. 
100-120 
channels is under design and construction at EMuS. 


At present, 
spectrometer with segmental 
The experimental studies of plastic scintillators, a 


single detection channel, a multi-channels 
spectrometer system and the data acquisition (DAQ) 


system are in progress. 
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